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Abstract 

A recently introduced nonlinear Fokker-Planck equation, derived directly from a master equa- 
tion, comes out as a very general tool to describe phenomenologically systems presenting complex 
behavior, like anomalous diffusion, in the presence of external forces. Such an equation is char- 
acterized by a nonlinear diffusion term that may present, in general, two distinct powers of the 
probability distribution. Herein, we calculate the stationary-state distributions of this equation 
in some special cases, and introduce associated classes of generalized entropies in order to satisfy 
the H-theorem. Within this approach, the parameters associated with the transition rates of the 
original master-equation are related to such generalized entropies, and are shown to obey some 
restrictions. Some particular cases are discussed. 
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I. INTRODUCTION 



The standard statistical-mechanics formalism, as proposed originally by Boltzmann and 
Gibbs (BG), is considered as one of the most successful theories of physics, and it has en- 
abled physicists to propose theoretical models in order to derive thermodynamical properties 
for real systems, by approaching the problem from the microscopic scale. Such a prescrip- 
tion has led to an adequate description of a large diversity of physical systems, essentially 
those represented by linear equations and characterized by short-range interactions and/or 
short-time memories. Although BG statistical mechanics is well formulated (under certain 



Drium systems, 
31. One of the 



restrictions) for systems at equilibrium, the same is not true for out-of-equi 
in such a way that most of this theory applies only near equilibrium |2 
most important phenomenological equations of nonequilibrium statistical mechanics is the 
linear Fokker-Planck equation (FPE), that rules the time evolution of the probability distri- 
bution associated with a given physical system, in the presence of an external force field {4], 
provided that the states of the system can be expressed by a continuum. This equation 
deals satisfactorily with many physical situations, e.g., those associated with normal diffu- 
sion, and is essentially associated with the BG formalism, in the sense that the Boltzmann 
distribution, which is usually obtained through the maximization of the BG entropy under 
certain constraints (the so-called MaxEnt principle), also appears as the stationary solution 



of the linear FPE 



Nevertheless, restrictions to the applicability of the BG statistical-mechanics formalism 
have been found in many systems, including for instance, those characterized by nonlinear- 
ities, long-range interactions and/or long-time memories, which may present several types 
of anomalous behavior, e.g., stationary states far from equilibrium [6,l7, f|. These anoma- 
lous behaviors suggest that a more general theory is required; as a consequence of that, 



many attempts have 
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j een made, essentially by proposing generalizations of the BG entropy 
13, Q]- Among these, the most successful proposal, so far, 
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appears to be the one suggested by Tsallis [9|], through the introduction of a generalized 
entropy, characterized by an index q, in such a way that the BG entropy is recovered in the 
limit q — > 1. The usual extensivity property of some thermodynamic quantities holds only 
for 9=1, and if q ^ 1 such quantities do not increase linearly with the size^of^he system; 
this has led to the so-called nonextensive statistical-mechanics formalism 
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Among many systems that present unusual behavior, one should mention those char- 
acterized by anomalous diffusion, e.g., particle transport in disordered media. A possible 
alternative for describing anomalous-transport processes consists in introducing modifica- 
tions in_the standard FPE. Within the most common procedure, one considers nonlinear 

;he cases come out as simple phenomenological generalizations 



FPEs 
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of the usual linear FPE 
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291 ] . In these nonlinear systems, 



interesting new aspects appear, leading to a wide range of open problems, in such a way that 
their investigation has led to a new research area in physics, with a lot of applications in 
natural systems. It is very common to find the power-law-like probabilit y d istribution that 
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maximizes Tsallis' entropy as solutions of some nonlinear FPEs 
seems that the nonextensive statistical mechanics formalism appears to be intimately related 
to nonlinear FPEs, motivating an investigation for a better understanding of possible con- 



nections between generalized entropies and nonlinear FPEs [19|, |20j, |21|, |25|, |29|, |31|, [32j, l33l. |34|. 

Recently, a general nonlinear FPE has been derived directly from a standard master 
equation, by introducing nonlinear contributions in the associated transition probabilities, 
leading to [35|, [36] 



dt dx dx V dx 



A) 



where a and b are constants, whereas /i and v are real exponents [38]. The system is in 
the presence of an external potential <j)(x), associated with a dimensionless force A(x) = 

X 

—d(j)(x) I dx [4>{x) = — f A(x')dx']; herein, we assume analyticity of the potential <p(x) and 

— oo 

integrability of the force A(x) in all space. 

In what concerns the functional Q[P(x, t)], we are assuming its differentiability and in- 
tegrability with respect to the probability distribution P(x,t), in such a way that at least 
its first derivative exists, i.e., that it should be at least fi[P] G C 1 . Furthermore, this func- 
tional should be a positive finite quantity, as expected for a proper diffusion-like term; this 
property will be verified later on, as a direct consequence of the H-theorem. 

As usual, we assume that the probability distribution, together with its first derivative, 
as well as the product A(x)P(x, t), should all be zero at infinity, 
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P(x,t)U ±oo = 



dP(x,t) 



dx 



A(x)P(x,t)\ x . 



• ±oo 



(Vt) . (2) 



x— >±oo 



The conditions above guarantee the preservation of the normalization for the probability 
distribution, i.e., if for a given time to one has that dx P(x,to) = 1, then a simple 
integration of Eq. ([1]) with respect to the variable x yields, 



| jH dz P(x, = - L4(*)P(x, t)}^ + f)] 



dP(x,t) 
dx 







(3) 



and so, 



/oo /»oo 
dx P(x,t) = / dxP(x,t ) = l (Vt) . 
-oo J — oo 



(4) 



In the present work we investigate further properties of the nonlinear FPE of Eq. (JTJ, finding 
stationary solutions in several particular cases, and discussing its associated entropies, that 
were introduced in order to satisfy the H-theorem. In the next section we present stationary 
solutions of this equation; in section III we prove the H-theorem by using Eq. (jTJ), and show 
that the validity of this theorem can be directly related to the definition of a general entropic 
form associated with this nonlinear FPE. In section IV we discuss particular cases of this 
general entropic form and their associated nonlinear FPEs. Finally, in section V we present 
our conclusions. 



II. STATIONARY STATE 



The nonlinear FPE of Eq. (JTJ) is very genera 



one related to Tsallis' thermostatistics 
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and covers several particular cases, e.g., the 



30|. In this section we will restrict ourselves 



to a stationary state, and will derive the corresponding solutions for particular values of the 
parameters associated with this equation. 

Let us then rewrite Eq. (Tj[|) in the form of a continuity equation, 



dP(x,t) dj{x,t) 



at 



dx 



0; j(x,t)=A{x)P(x,t)-n[P(x,t)] 



dP(x,t) 
dx 



(5) 
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in such a way that a stationary solution of Eq. (pQ), P st (x), is associated with a stationary 
probability flux, j s t{x) = constant, which becomes j s t{x) = 0, when one uses Eq. ([2]). 
Therefore, using the functional £l[P] of Eq. flTJ, the stationary-state solution satisfies, 



A(x) = [a^-\x) + 6(2 - u)P^ 2 (x)] 



dP st {x) 
dx 



(6) 



which, after integration, becomes 



^-<P(x) = a-^-P^ 1 (x) 
H — 1 



(7) 



where 0o represents a constant. The equation above may be solved easily in some particular 
cases, e.g., v — ji, v — 2, and /i — 0, 
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(8) 



where a = /j, in the cases v = [i and z/ = 2, whereas a = i/, if /x = 0. In the equation above, 
M + — 2/) f° r V > 0, and zero otherwise. Another type of solution applies for /j = 2z/ — 1, or 
i/ = 2/i- 1, 



i 

a-1 



= / dx 



where 



i 

ct-l 



(9) 



(10) 



2b (3-2/x)(//-i; 



k a 



if i/ = 2/x — 1 (a = /i) 



and we are assuming that [1 + K(<f)(x) — </>o)] > 0. Some well-known particular cases of 
the stationary solutions presented above come out easily, e.g., from Eq. (jSj) one obtains, in 



all three situations that yielded this equation, the exponential solution associated with the 
linear FPE in the limit a — > 1 [with <fro oc (a — as well as the generalized exponential 

solution related to Tsallis thermostatistics, for a = 2 — q, where q denotes Tsallis' entropic 
index. 

III. H THEOREM AND THE ASSOCIATED ENTROPY 

In this section we will demonstrate the H-theorem by making use of Eq. (CD), and for 
that purpose, an entropic form related to this equation will be introduced. Let us therefore 
suppose a general entropic form satisfying the following conditions, 

oo 

/j2 
dx g[P{x)} ■ g[0] = ; g[l] = ; ^| < , (10) 

— oo 

where one should have g[P(x,t)} at least as g[P(x, t)] 6 C 2 ; in addition to that, let us also 
define the free-energy functional, 

oo 

F = U —~ S ; U = J dx <f>(x)P{x, t) , (11) 

— oo 

where /3 represents a Lagrange multiplier, restricted to (3 > 0. Furthermore, we will show 
that this free-energy functional is bounded from below; this condition, together with the 
H-theorem [(dF/dt) < 0), leads, after a long time, the system towards a stationary state. 

A. H Theorem 

The H-theorem for a system that exchanges energy with its surrounding, herein repre- 
sented by the potential <f)(x), corresponds to a well-defined sign for the time derivative of 
the free-energy functional defined in Eq. ( flTT) . Using the definitions above, 
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oo oo 
-oo — oo 



(12) 



Now, one may use the FPE of Eq. ([I]) for the time derivative of the probability distribution; 
carrying out an integration by parts, and using the conditions of Eq. ([21), one obtains, 



OF 
~dt 



-TT- — — I dx 



d(f>(x) ld 2 g[P]dP 
dx (3 dP 2 dx 



(13) 



Usually, one is interested in verifying the H-theorem from a well-defined FPE, together 
with a particular entropic form, in such a way that the quantities Q[P] and d 2 g[P}/dP 2 
are previously defined (see, e.g., Refs. JJJ, |37|]). Herein, we follow a different approach, 
by assuming that the general Eqs. ([1]) and (flOl) should be satisfied; then, we impose the 
condition, 



d 2 g[P} 

dP 2 



-p 



P{x,t) 



(14) 



in such way that, 



dF 
~dt 



dx P(x, t) 



d<p(x) Q[P] dP~ 
dx P(x, t) dx 



< 



(15) 



It should be noticed that Eq. (JHj) , introduced in such a way to provide a well-defined sign 
for the time derivative of the free-energy functional, yields two important conditions, as 
described next. 

(i) £l[P] > [cf. Eq. (II Op ], which is expected for an appropriate diffusion-like term. 

(ii) It expresses a relation involving the FPE of Eq. ([T|) and an associated entropic form, 
allowing for the calculation of such an entropic form, given the FPE, and vice-versa. Since 
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the FPE is a phenomenological equation that specifies the dynamical evolution associated 
with a given physical system, Eq. (I14p may be useful in the identification of the entropic 
form associated with such a system. In particular, the present approach makes it possible 
to identify entropic forms associated with some anomalous systems, exhibiting unusual be- 
havior, that are appropriately described by nonlinear FPEs, like the one of Eq. (1T1). As an 
illustration of this point, let us consider the simple case of a linear FPE, that describes the 
dynamical evolution of many physical systems, essentially those characterized by normal 
diffusion. This equation may be recovered from Eq. ([1]) by choosing // = v — 1, in such a 
way that £l[P] = a + b = D, where D represents a positive constant diffusion coefficient with 
units (time) -1 . One may now set the Lagrange multiplier (3 = ks/D, where ks represents 
the Boltzmann constant; integrating Eq. (JHJ), and using the conditions of Eq. (JTUJ), one gets 
the well-known BG entropic form, 

g[P} = -k B P\nP . (16) 

In the next section we will explore further the relation of Eq. (1141) . by analyzing other 
particular cases. 

The simplest situation for which condition (i) above is satisfied may be obtained by 
imposing both terms of the functional £l[P] to be positive, which leads to 

a > ; // > , or a < ; fj, < , (17) 

and 

b>0; v < 2 , or 6 < ; v > 2 . (18) 

It should be stressed that it is possible to have £l[P] > with less restrictive ranges for 
the parameters above. However, an additional property for the free-energy functional of 
Eq. (1111) to be discussed next, namely, the boundness from below, requires the conditions of 
Eqs. ( ITT)) and (TTH)) . with the additional restriction v > 0. 

Now, integrating Eq. (fTT)) for the general functional £l[P] of Eq. (JTJ, and using the 
standard conditions for g[P] defined in Eq. (TTOj) . one gets that 
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g[P] = -P 



-P" + b 



V 



v{y 



au(l -u) + 6(2 - - //) 
(l-fJ.)(l-u)u 



p 



(19) 



This entropic form recovers, as particular cases, the BG entropy (e.g., when ^t, v — > 1) and 
several generalized entropies defined previously in the literature, like those introduced by 

171 ] . Such particular cases, as 



n 



Tsallis |9j, Abe 10J, Borges-Roditi 12j, and Kaniadakis 16, 
well as their associated FPEs, will be discussed in the next section. 

For the simpler situation of an isolated system, i.e., 4>(x) = constant, the H-theorem 
should be expressed in terms of the time derivative of the entropy, in such a way that 
Eq. (fT3l) should be replaced by 



dS[P] 



dP\ (d 2 g[P]dP 



dP 2 dx 



dx Q[P] 



d 2 g[P] 
dP 2 



ox 



(20) 



In this case all that one needs is the standard condition associated with the FPE [same 
condition (i) above], i.e., > 0, and the general restrictions of Eq. (11 01) for the entropy. 
A similar result may also be obtained by proving the H-theorem using the master equation 
from which Eq. ([T]) was derived, with the transition probabilities introduced in Ref. 
(see the Appendix). 



35 



36] 



B. Boundness from Below 

Above, we have proven that the free-energy functional decreases in time, and so, for the 
existence of a stationary state at long times of an evolution process, characterized by a 
probability distribution P st (x), one should have that 



F(P(x,t))>F(P st (x)) (Vt). 



(21) 



In what follows, we will show this inequality and find the conditions for its validity. There- 
fore, using Eqs. ([7]), (fTUI) . and (TTTj) . we can write, 
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oo oo 

F(P(x,t))= J P(x,t) (t - a ^P!T\x)-b^P:f 1 (xfjdx-± J g[P]dx, (22) 



and so, 



F(P st )-F(P)= J (P-P st ) (a-J^ + b^P^ dx + \ j (g[P]-9[Pst})dx, 

— CO — oo 

(23) 

where we have used the normalization condition for the probabilities. Now, we insert the 
entropic form of Eq. (lT§j) in the equation above to obtain, 



oo 

F(P st )-F(P)= J [^Ps^.iP/Pstl + ^^P^iP/Pst] 



dx , (24) 



where, 

T a [z] = 1 — a + az — z a (a = /i, u) . (25) 

By analyzing the extrema of the functional T a [P/P st ], one can see that T a [P/P st ] > for 
< a < 1 and T a [P/P st ] < for a > 1 and a < 0. Therefore, the inequality of Eq. (12~TT) 
is satisfied for the range of parameters specified by Eqs. ( 1171) and ( 1181) . if one considers 
additionally, v > 0. 

For the case of an isolated system, the stationary solution turns out to be the equilibrium 
state, which is the one that maximizes the entropy. Therefore, one may use the concavity 
property of the entropy [cf. Eq. (TTTJ1) ] in order to get, 

oo 

S[P eq (x)]-S[P(x,t)] = J (g[P eq ]-g[P])dx>0 (Vt). (26) 
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IV. SOME PARTICULAR CASES 



In this section we analyze some particular cases of the entropic form of Eq. ([19]) and using 
Eq. (1141) . we find for each of them, the corresponding functional Q[P] of the associated FPE. 
In the examples that follow, we will set the Lagrange multiplier j3 = k/D, with k and D 
representing, respectively, a constant with dimensions of entropy and a constant diffusion 
coefficient. 

(a) Tsallis entropy 9]: This represents the most well-known generalization of the BG entropy, 
which has led to the development of the area of nonextensive statistical mechanics 6|, |7|, |8j . 
One may find easily that Eq. fTITJT) recovers Tsallis entropy in several particular cases, e.g., 
{b = 0, a = D,fi = q}, {a = 0,b = Dv/(2-v),v = q}, and {a = D/2,b = Du/[2(2-u)], /i = 
v = q}. For all these cases one may use Eq. (fl4j) . in order to get the corresponding functional 
Q[P], 

g[P] = k^^, n[P] = qDP*- 1 . (27) 



With the functional £l[P] above, one identifies the nonlinear FPE that presents the well- 
known q -exp onential, or Tsallis distribution (replacing q 2 — q), as a time-dependent 



solution 120 
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25] 



(b) Abe entropy 10j: This proposal was inspired in the area of quantum groups, where 
certain quantities, usually called g-deformed quantities, are submitted to deformations and 
are often required to possess the invariance q <-> g" 1 . The Abe entropy may be obtained from 
Eq. ffl~9"l) in the particular case {a = D(q — l)/(g — q^ 1 ), b = —Dq(q + l)/[(q — q~ 1 )(q + 2)},^ = 
— v = g}, for which 



p? _ p-g 

g[P] = -k r , Q[P] 

q-q 



D 



g(g~ jO pq-l _ g(g+ 1 ) p- (? -l 



q-q- 



q-q' 



(28) 



(c) Borges-Roditi entropy [12|: This consists in another generalization of Tsallis entropy, 
where now one has two distinct entropic indices, q and g', with a more general invariance 
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q <-> q'] this case may be obtained from Eq. (fT9l) by choosing {a = D(q — l)/(q' — q),b 
Dq'{q> - l)/[{q - q')(2 - q')],fi = q,v = q'}. One gets, 



g[P\ 



-k 



pi _ pi 



q-q' 



* , 0[P] = D-!—; (q(q - l)P^ - q'(q> - 1)P^ 



(29) 



(d) Kaniadakis entropy 16|, ll7(|: This is also a two-exponent entropic form, but slightly 
different from those presented in examples (b) and (c) above; it may be reproduced from 
Eq. ffT9"j) by choosing {a = b = D/[2(l + q)], fi = 1 + q, v — 1 — q}, in such a way that 



g\P] = -—( -^—P 1+q l —P l - q \ , Q\P] = -(P q + P~ 

y[ 1 2q\l + q 1-q J 2 



(30) 



Except for the well-known example (a), i.e., Tsallis entropy and its corresponding FPE, 
the other three particular cases presented herein were much less explored in the literature. 
Their associated FPEs, defined in terms of their respective functionals £i[P] above are, 
to our knowledge, presented herein for the first time. These equations, whose nonlinear 
terms depend essentially in two different powers of the probability distribution, may be 
appropriated for anomalous- diffusion phenomena where a crossover between two different 



diffusion regimes occurs 



24] 



V. CONCLUSION 



We have analyzed important aspects associated with a recently introduced nonlinear 
Fokker-Planck equation, that was derived directly from a master equation by setting non- 
linear effects on its transition rates. Such equation is characterized by a nonlinear diffusion 
term that may present two distinct powers of the probability distribution; for this reason, 
it may reproduce, as particular large range of nonlinear FPEs of the literature. 

We have obtained stationary solutions for this equation in several cases, and some of them 
recover the well-known Tsallis distribution. We have proven the H-theorem, and for that, 
an important relation involving the parameters of the FPE and an entropic form was intro- 
duced. Since the FPE is a phenomenological equation that specifies the dynamical evolution 
associated with a given physical system, such a relation may be useful for identifying en- 
tropic forms associated with real systems and, in particular, anomalous systems that exhibit 
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unusual behavior and are appropriately described by nonlinear FPEs. It is shown that, in 
the simple case of a linear FPE, the Boltzmann-Gibbs entropy comes out straightforwardly 
from this relation. Considering the above-mentioned nonlinear diffusion term, this relation 
yields a very general entropic form, which similarly to its corresponding FPE, depends on 
two distinct powers of the probability distribution. Apart from Tsallis' entropy, other en- 
tropic forms introduced in the literature are recovered as particular cases of the present one, 
essentially those characterized by two entropic indices. Nonlinear FPEs (as well as their 
associated entropic forms) whose nonlinear terms depend essentially on two different powers 
of the probability distribution, like the ones discussed in the present paper, are good candi- 
dates for describing anomalous- diffusion phenomena where a crossover between two different 
diffusion regimes may take place. As a typical example, one could have a particle trans- 
port in a system composed by two types of disordered media, characterized respectively, by 
significantly different grains (e.g., different average sizes), arranged in such a way that the 
diffusion process is dominated by one medium, at high densities, and by the other one, at 
low densities. 
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Appendix 

In this appendix we will prove the H-theorem directly from the master equation, for an 
isolated system (i.e., no external forces). Let us consider a system described in terms of W 
discrete stochastic variables; then, P n (t) represents the probability of finding this system in 
a state characterized by the variable n at time t. This probability evolves in time according 
to a master equation, 

dP n W 

7^ = ^ ^ {Pm w m,n Pn w n,rri) > (31) 
m=l 

where w k ,i(t) represents the transition probability rate from state k to state I. The nonlinear 



effects were introduced through the transition rates 35| 
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w k>l (A) = i,(S k>l+1 + 5 k ^ x )[aPr\t) + bPr\t)} , (32) 



where A represents the size of the step of the random walk. Herein, we shall consider a 
random walk characterized by A = 1; substituting such a transition rate in the master 
equation one gets, 

BP 

n = 1 : -^± = aP£- aP? + bP 2 P^ 1 - bP x P^ x , (33) 
BP 

n = W: -^- = aP^_ x - aP^ + bPw^P^ 1 - bP w P^\ , (34) 
BP 

n = 2, (W — 1) : = a (P% +1 + - 2aP£ + bP^ 1 (P n+1 + P n _0 

-bP n (P^l + P»:l) , (35) 



where we have treated the borders of the spectrum (n — 1 and n = W) separately from the 

w 

rest. Let us now consider the entropy, S = ^ <7[-fn]> satisfying the same conditions specified 

71=1 
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in the text [see Eq. (fIDI) ]; the H-theorem, to be proven below, is expressed in terms of its 
time derivative, 



w 



dS i^dg[P n ]dP n n 
— = > 7 J / > . 36 
dt ^ dP n dt ~ y ' 

n=l 



Then, using Eqs. fl33l) - (l35|) one can write this time derivative as, 



w-i 



n=l 
W 



+ E ^ ( aP n~l ~ aP X + bPn-lPr 1 ~ bPnP^l) 



(37) 



n=2 



W-l 

The sum indices can be rearranged to yield all summations in the range 5 an d thus the 

n=l 

time derivative of the entropy becomes, 



ds 

dt 



E - 1^ 1 ) [« - Pit) - bp n p n+ , (k-j - pr)] • 

n=l ^ n n+1 J 



(38) 



The negative curvature of the entropic function [cf. Eq. (flOl) ] implies that its first derivative 
decays monotonically with P n . Hence, for P„+i > P n , the condition of Eq. fl36l) is satisfied 
for 

pM _ p/i pv-2 pi/— 2 

a p n+1 " - bP n P n+1 ^±± > , (39) 

-< n+1 — Pi n+1 — " n 

where we divided the term inside the brackets in Eq. (l3"8l by the difference AP = P n+ \ — P n . 
Now we consider the limit AP — > and obtain, 

a/iP^ 1 + 6(2 - v)P v ~ 1 > , (40) 



which corresponds to the condition fi[P] > found in the text, when proving the H-theorem 
by making use of the FPE of Eq. ([1]). It should be mentioned that the procedure above works 

1(3 



also for P n > P n +i ', therefore, in this appendix we have proven that the H-theorem holds 
for an arbitrary state n and all times t of an isolated system. 
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